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ABSTRACT 

We present a detailed analysis of high resolution near-infrared imaging and spectroscopy of the 
potential star cluster IRS13E very close to the massive black hole in the Galactic Center. We detect 
19 objects in IRS13E from Ks-band images, 15 of which are also detected reliably in H-band. We derive 
consistent proper motions for these objects from the two bands. Most objects share a similar westward 
proper motion. We characterize the objects using spectroscopy (1.45 to 2.45 /ini) and (narrow-band) 
imaging from H- (1.66 ^m) to L'-band (3.80 /im). Nine of the objects detected in both Ks- and H-band 
are very red, and we find that they are all consistent with being warm dust clumps. The dust emission 
may be caused by the colliding winds of the two Wolf-Rayet stars in the cluster. Three of the six 
detected stars do not share the motion or spectral properties of the three bright stars. This leaves 
only the three bright, early- type stars as potential cluster members. It is unlikely that these stars 
are a chance configuration. Assuming the presence of an IMBH, a mass of about 14000 7\f0 follows 
from the velocities and positions of these three stars. However, our acceleration limits make such an 
IMBH nearly as unlikely as a chance occurrence of such a star association. Furthermore, there is no 
variable X-ray source in IRS13E despite the high density of dust and gas. Therefore, we conclude 
that is unlikely that IRS13E hosts a black hole massive enough to bind the three stars. 
Subject headings: infrared: stars — Galaxy: center 



1. INTRODUCTION 

The innermost parsec of the Galaxy hosts the su- 
permassive black hole (SMBH) Sgr A* (jSchodel et al.l 
120021 : iGhez et all l2003| ) acc ompanie d by a population 
of young WR/O-stars (HbTrest et all 119871: lAllen et al.l 
\T9m IKrabbe et all I199H IGenzel et al.l \199d\ . Most 
of these stars reside in one or two disk-like structures 
(IGenzel et al.ll2003l: iPaumard et al.|[200llLu et al|[200l 
iBartko et al.ll2009[ ). A group of at least three such bright 
stars called IRS13E at a distance of 3.5" = 0.13 pc from 
Sgr A* is of special in terest. This gro up has a diame- 
ter of about 0.5". Maillard et al.l (j2004D identified two of 
the stars, E2 and E4, from their emission lines as early- 
type stars. In addition they identified four other early- 
type stars in IRS13E from broad-band SED fitting. The 
four brightest stars share a c ommon proper motion (|Ottl 
I2002D . .Maillard et al.l (|2004D concluded without using a 
statistical test that such an association of young stars 
cannot be a coincidence. From the radial velocities of 
two of the stars they estimated a cluster mass of at least 
750 AIq if the two stars are bound. This mass is higher 
than the stellar mass seen in IRS13E. They explained the 
additional mass by the presence of an intermediate mass 
black hole (IMBH). 

In the simulations o f ste llar clusters of 
iPortegies Zwart fc McMillanI ()2002[ ) the core of a 
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dense cluster collapses and forms an IMBH. Such an 
IMBH would be necessary, if IRS13E were an inspiraling 
cluster that survives the infall into the Galactic Center 
(GC ) and reaches the central pars ec before disintegrat- 
ing |Hansen &: Milosavlievid [2003D . A cluster without 



33|). 

central mass would be disrupted (|Gerhardll2001f ) by the 
ti dal forces of the SMB H. 

iSchodel et al.l (|2005D measured the proper motions of 
the four brightest sources more accurately and estimated 
that the cluster has a mass of about 10000 to 50000 Mq if 
it is gravitationally bound. According to these authors, 
an IMBH of this mass is unlikely, mainly because of the 
lack of radio and X-ray emission in IRS13E, despite the 
presence of a lot of dust. They suggested that IRS13E 
could also be a cluster in the process of dissolution or a 
c hance association . 

IPaumard et al.l ()2006[ ) identified the spectral types of 
the three brightest early-type stars. In addition, they 
measured the stellar surface density around IRS13E on a 
deconvolved H-band image. Inside a radius of 0.30", the 
core of IRS13E, they found at least twelve stars. Further- 
more, they found an over-density out to 0.68" which has a 
total significance of 4.5 a. They concluded that IRS13E 
is a cluster. Given that velocities had been m easured 
for only four of the stars, IPaumard et al.l (|2006[ ) argued 
that the total velocity dispersion could be small. Then 
no dar k mass is needed to explain IRS13E. iTrippe et al.l 
(|2008[ ) measured the velocities of the stars between 0.30" 
and 0.68". The velocities of most of these stars are dif- 
ferent from the stars in the center. Therefore most of 
them are not cluster members and thus the overden sity 
is less significant compared to IPaumard et al.l (|2006l ) . 

IRS13E2 and IRS13E4 are ca ndidate members of 
the face-on counte rclockwise disk ([Paumard et al.l [20061 : 
IBartko et al.l l2009l) . In this case the physical distance 
of these stars to Sgr A* is identical with the projected 
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one. According to iPaumard et"all (|200l [20061 ) IRS13E 
is embedded in the b ar of the mi nispiral at a distance of 
z= 7 - 20" of Sgr A* (|Lisztil2003D . In this case IRS13E is 
rather far away from Sgr A* and IRS13E2 and IRS13E4 
would not be candidates members of the counterclock- 
wise disc (Bartko et al. 2009 ). 

Accordingly, the nature of IRS13E is still a matter of 
debate. This is largely because the proper motions and 
the nature of the member objects are known only for very 
few sources in the core with a radius of about 0.3". In 
particular, it is not clear whether all objects identified 
by previous works actually are stars. In this paper we 
analyze seven years of imaging data (presented in sec- 
tion 2). We identify 19 objects in the core of IRS13E, 
for all of which we can measure proper motions (section 
3). For the brightest sources El and E2, we obtain ac- 
celeration limits which start to to constrain the nature 
of IRS13E (section 3). From H- to L'-band photometry 
and H-|-K-band spectra we constrain the nature of the 
objects (section 4). In Section 5 we calculate the proba- 
bility that our acceleration limits occur for the two cases 
that a) IRS13E is bound by an IMBH and b) that it is 
a chance association. Section 6 considers other data in a 
more qualitative fashion and theoretical models. Finally 
we summarize and conclude in Secti on 7. We assume 
a distance to the GC of Rq = 8kpc (|Reidl fl993l) and a 
mass of the SMBH of Af = 4 x lO^M© (|Ghez et al.ll200l 
iGillessen et al...2009 ). 

2. DATASET 

This study uses seven years of adaptive optics based 
imaging data and spectroscopic data obtained with an 
integral field spectrograph in 2004 and 2009. This section 
briefly describes our data set. 

2.1. NACO 

We obtained images with NAOS /CONICA (NACO) 
mounted on UT4 at the VLT (Lenzen et all l2003t 
iRousset et all [200I FI. The camera CONICA, together 
with the adaptive optics (AO) system NAOS, achieves 
diffraction limited images in the near-infrared. For this 
study, we use mainly the 13 mas/pixel image scale. On 
this scale, the point spread function (PSF) is oversam- 
pled. This is suited better for highly crowded fields than 
a strict Nyquist sampling (jTrippe et al.ll2009| ). 

For source detection and astrometry, we use data from 
nearly every observing run from April 2002 to May 2009, 
excluding only the worst 10% of the observing runs by 
eye. In total, our sample consists of 37 H-, 74 Ks- and 5 
L'-band images. 

We base the broad-band photometry on four H-, four 
Ks- and three L'-band images that were of particularly 
good quality. Furthermore, we use four narrow-band 
images for photometry, which we obtained at 2.06 /im, 
2.24 ^m (12-06-2004), 2.17 fim and 2.33 ^im (13-06-2004). 

All images are flat-fielded, bad pixel corrected and sky 
subtracted. In addition, we correct images in which sat- 
uration is import ant for the nonlinearity of the detector 
()Fritz et al.ll2010D . 

^ We use NACO data from ESO programs 70.B-0649, 71.B-0077, 
71.B-0078, 71.B-0365, 072.B-0285, 073.B-0084, 073.B-0085, 073.B- 
0665, 073.B-0775, 075.B-0093, 077.B-0014, 078.B-0136, 179.B- 
0261, 179.B-0932, 183.B-0100. 




Fig. 1. — Lucy-deconvolved Ks-image of the GC from 13-03-2008, 
smoothed with a Gaussian with a FWHM of 4 pixels. The intensity 
scale is logarithmic. Sgr A* and IRS13E are marked. 



We combine single images of typically 30 seconds ex- 
posure into stacked images. Because the faint additional 
sources in IRS13E lie in the PSF halos of the bright stars, 
the prime goal of the data reduction is to obtain the high- 
est Strehl ratio possible fo r each data set, e ven at the cost 
of smaller total exposure (jFritz et al.ll2010l) . For this rea- 
son, we select frames of good quality (small FWHM) for 
a combined image of an observing run only. 

2.2. SINFONI 

Spectroscopy allows one to obtain a secure identi- 
fication of spectral types and to measure radial ve- 
locities. We use spectra obtained wit h the integral 
field spectrom eter SINFONI (Eisenhau er et al.l 120031 : 
iBonnet et all [20031 at UT4 of the VLTQ. We analyze 
high-quality H-|-K-band laser and natural guide star data 
from May 21 and 23, 2009. The data from the latter 
date have a Strehl ratio of 28% at 2.11 /xm and a total 
of 2400 s exposure. The spatial sampling is 12.5 mas x 
25 mas pixel"^, the spectral resolution is 1500. We apply 
the standard data reduction for SINFONI data, includ- 
ing detector calibrations (such as bad pixel correction, 
flat-fielding and distortion correction) and cube recon- 
struction. The wavelength scale is calibrated by means 
of emission line lamps and flnetuned with atmospheric 
OH lines. The remaining uncertainty corresponds to typ- 
ically (less than) lOkm/s. 

3. ASTROMETRY 
3.1. Image processing 

One of the main aims of this study is to characterize the 
faint objects in the core of IRS13E, the velocities of which 
are unknown. We base our analysis on deconvolution of 
the images, which requires the determination of the PSF 
for each image. 

For creating a P SF we use starfinder version 1.2 
piolaiti et al.ll2000D . We extract the PSF from ten to 15 

^ We use SINFONI data from ESO program 183.B-100. 
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stars brighter than mxs = 13.5 and symmetrically dis- 
tributed around IRS13E up to a dista nce of 4". W e then 
use the Lucy- Richardson algorithm (jLucvl 119741) . De- 
pending on the image quality, we choose between 6000 
and 20000 iterations. We are confident that the LR al- 
gorithm applied in this way also mostly correctly han- 
dles ext ended sources. Fo r example, the extended object 
IRS2L (Hi um et all |1996[) contains only one maximum 
after deconvolution, and a large fraction of the flux re- 
mains detached from the central maximum. Relatively 
faint sources which are in the same distance from a bright 
star as the radius of the PSF are sometimes influenced 
by deconvolution rings. Most (and all important) sources 
are either clearly outside or clearly inside of these rings. 
We do not smooth the result of the deconvolution at all 
and find stellar positions by centroids around local max- 
ima. 

We deconvolve all H- and Ks-band images. We then 
select by visual inspection the images for the subsequent 
analysis, on which at least some of the fainter objects 
are visible. These are 61 Ks-band images and 18 H- 
band images. This subjective quality selection is nearly 
identical to a Strehl ratio cut of 9% in both bands. 

3.2. Object detection 

Some care is needed to distinguish between real and 
spurious sources since we are dealing with a crowded 
field imaged at the diffraction limit. Our procedure is 
as follows. We start by measuring the positions of the 
local maxima on the quality-selected, deconvolved im- 
ages. We treat each frame individually, such that we can 
neglect maxima that are unlikely to correspond to a real 
source. These are for example maxima only marginally 
brighter than the noise floor or the background from the 
minispiral, or maxima that apparently are affected by 
deconvolution artifacts. 

We consider an object to be well-detected if it is found 
in either at least three images in both H- and Ks-band, or 
in at least ten Ks-band images. This is a robust criterion, 
since the objects are present in images that sample a 
wide range of PSF shapes. According to these criteria, 
we detect 19 sources in the core (r < 0.3") of IRS13E 
(figure [5]). The rms position deviation of these sources 
to a linear fit after transformation into the astrometric 
reference frame to the Ks-band data is 7 mas in median. 
Thus all measured positions of the 19 sources lie within 
an area of 4000 mas^, only 1 /50 of the area of the cluster 
core. Therefore the sources are not random fiuctuations. 

We test by means of simulated images that indeed our 
detection procedure is robust. Using a variety of PSFs 
we find that it is only possible to continue to obtain the 
same sources after deconvolution if indeed these sources 
had been put into the simulation. 

3.3. Proper motions 

Given the pixel positions, we go to an astrometric co- 
ordinate system by constructing a full linear transforma- 
tion u sing ten brigh t , isola ted stars from the astrometric 
list of iTrippe etall (|2q08[ ) close to IRS13E. Proper mo- 
tions are in turn obtained by fitting the positions as a 
function of time by linear functions. Since we are using 
many position measurements for the fitting, we simply 
use an unweighted fit. The errors for any given object are 



Fig. 2. — Coadd of four Lucy-deconvolved Ks-images of IRS13E 
between 2002 and 2007, smoothed with a Gaussian with a FWHM 
of 2 pixels. The intensity scale is logarithmic. The sources which 
wc consider real are marked with numbers. Nearly all of these 
objects are also found in H-band images. 

set to the scatter of the data around the linear motion 
model. For the fit, 5tT-outliers are rejected. There are 
only two such cases in H and Ks for all the objects. The 
absolute position towards Sgr A* has an additional un- 
certainty of about 2 mas CR eid et al.ll2007t iTrippe et al.l 
[2008; Gillcssen et al. 2009) corresponding to a relative 
error of « 0.05% which safely can be neglected in our 
analysis. 

In the L'-band the exclusion of 5(T outliers becomes 
important, and the errors are mostly larger. The L'-band 
velocities differ significantly from the velocities obtained 
from the other bands for the three brightest stars. The 
differences can be as large as 16(t and 2.5mas/yr. The 
differences are probably caused by L'-band dust emission 
that is important around IRS13E. The L'-band velocities 
are less useful and we neglect them in our analysis. 

The proper motion fits are presented in table [T] and 
figure [T51 For the bright stars El, E2 and E4.0 we also 
fit second order polynomials to the stellar positions as a 
function of time to test for accelerations. The accuracy 
reached for the three brightest stars starts to constrain 
the mass of a putative IMBH. 

For most objects the velocities derived from H- and 
Ks-band data are consistent. The biggest discrepancy 
happens for E4.0 that shows a velocity difference of 16 
km/s or 3.7 cr in R.A. The value for the R. A. -velocity in 
Ks lies between the velocities obtained from H- and L'- 
band data. Therefore, we assume that dust perturbs the 
position of this star also in Ks-band, not only in L'-band. 
Consequently we only use the astrometry from H-band 
data for this star, and we also don't trust the marginally 
significant acceleration in Dec. direction. 

Similarly, for the faint or confused objects (E3, E4.1, 
E4.2, E5, E7, E8) we adopt the measurements derived 
from the band where the source is more prominent com- 
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pared to the local background level. For the relatively 
blue objects E3.5, E4.0 and E4.2 this is the H-band, for 
the other ones Ks-band. For the unconfused objects El, 
E2 and E6 we use the weighted average of the two bands. 

4. CHARACTERIZATION OF THE OBJECTS 

In order to address the question of whether IRS13E 
hosts an IMBH, it is important to know the nature of 
its constituent objects. Only stars can be used as fidu- 
cial tracers for a gravitational potential, and an over- 
density of objects compared to the stellar background 
is only meaningful if these are actually stars. In this 
section we first describe the combined photometric and 
spectroscopic method to construct a SED. Then we use 
it to constrain the nature of the individual objects. We 
present the results for the most important objects here 
and leave for the appendix a description of the details of 
the method and the results for the remaining objects. 

4.1. Method 

For each object wc construct a spectral energy distri- 
bution (SED), ideally ranging from H- to L'-band. For 
this we use broad- and narrow-band images plus the two 
spectro-imaging data cubes in the combined H-f K-band. 

Where possible, we identify the objects in the SIN- 
FONI data by their spectral features as late- or early- 
type stars. In this way we show that El, E2, E4.0 and 
E6.0 are stars. For the other objects we construct SEDs 
from a deconvolved SINFONI cube and the narrow-band 
imaging data (appendix |A|). We then estimate whether 
a stellar component is requir ed to explain th e SEP . 

Following the approach of iMaillard et al.l (|2004D , we 
consider models composed of an extinction parameter 
and either one or two temperatures corresponding to a 
star and dust. We derive the temperature of the stellar 
component from the lines in the spectrum when possible. 
Otherwise we assume a temperature of 9480 K (corre- 
sponding to the spectral type AOV) for the star. This is 
sufficient, since the shape of the SEDs in our wavelength 
range depends only weakly on the assumed temperature 
for the stellar component of the SED model. 

We use the extinction law of .Drain^ (119891) sho rtward 
of Brackett-/? and Aks/Al' = 1.75 |Lutz et al.llT996l ). 
When fitting models to the less well measured SEDs, we 
only use extinction values that differ less than AAks = 
0.3 from the value for El. This range seems adequate 
given the extinction map of Schodel et al. (2010) and our 
result does not depend strongly on this assumption. We 
neglect the nebular lines originating from the interstellar 
gas. At our precision their contributions to the flux are 
negligible 

For the bright early-type stars we have well-measured 
spectra such that an explicit comparison with atmo- 
sphere models becomes viab le. We use the code CM- 
FGEN (jHilher fc Milleii[T998l) to derive the main param- 
eters of the stellar atmospheres of El, E2 and E4.0. 
The code computes non-LTE atmosphere models in- 
cluding wind and li ne-blanketing. We proceed as in 
iMartins et al.l (|2007[ ). A slight improvement in the 
method is the use of iterations on the hydrodynamic 
structure to ensure a better consistency between the at- 
mosp here model per se a nd the density /velocity struc- 
ture (jMartins et al.l [2(309 1) . The derived stellar parame- 
ters rely on line ratios and photometry. In addition, we 
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Fig. 3. — H-l-K-band spectrum and photometry of the two WR- 
stars in IRS13E: E2 (WN8) and E4.0 (WC9). We fit full atmo- 
sphere models plus blackbodies with dust-like temperatures to the 
observed spectra and L'-band fiux points. The asterisks denote the 
average flux of the models in L'-band, to which the measured flux 
needs to be compared. 

derive the extinction by fitting the shape of the infrared 
SEDs. Hence, for each star, we have computed several 
models with different luminosities and looked for the best 
combination of luminosity, extinction and dust emission. 

4.2. Individual Objects 

This section describes the nature of the most important 
objects in IRS13E. For a description of the other objects, 
see appendix [B1 Table [5] presents the values of the model 
parameters for all objects, for which our data allow an 
identification. 

4.2.1. IRS13E3.0 

We resolve E3 on the deconvolved images into six sta- 
ble sources and name them E3.0, E3.1, E3.2, E3.3, E3.4 
and E3.5. This section describes the brightest compo- 
nent E3.0, which likely is extended or multiple at our 
resolution. Its flux increases signiflcantly by 0.25 mag 
between 2002 and 2008 in H- and Ks-bands, the main 
increase happens between 2004 and 2006. This is not 
caused by confusion. 

We extract a spectrum from the SINFONI data at the 
position of E3.0 (flgureH]). This spectrum also contains 
flux from the other five E3 sources and the extended 
background, but E3.0 completely dominates due to its 
brightness. 

No stellar line (absorption line, or broader emission 
line) is visible in E3.0, but many narrow emission lines 
from H, Fe and He are present. Nearly all these 
lines are very similar to the surrounding nebular back- 
ground. A few lines like H2 at 2.122 /xm and 2.223 /im 
are only present in the background. Most identifiable 
line s (HI, Hel, Fell, Felll) are known from the minispi- 
ral (|Lutz et al.lll993l ) or are other hydrogen lines. 

Stellar absorption lines of fainter stars are likely too 
weak for being detectable against the strong narrow emis- 
sion lines of the surrounding gas. However, that there are 
any WR-stars in E3.0 seems unlikely, since the strong, 
broad emission lines of such stars would be visible. For 
example, E4.0 (a WC9 star) scaled down to toks = 15 
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TABLE 1 

Proper motions and acceleration limits of the objects in IRS13E 



Object 


band 


R.A. 


VR.A. 


Dec. 


VDcc. 


flR.A. 


aDcc. 






[as] 


[km / s] 


[as] 


[km / s] 


[^as/yr^] 


[^as/yr^] 


El 


ti and Ks 


-2.958 


-142.5 ± l.o 


-1 


645 


-1U5.5 ± U.7 


22 it 20 


11 it 14 


E2 


ti and Ks 


-3.171 


-249.2 ih 1.1 


-1 


732 


23 it 0.7 


20 it 19 


-14 it 16 


E3.0 


Ks 


-3.184 


-82 ± 9 


-1 


528 


8 it 9 






E3.1 


Ks 


-3.126 


19 ± 39 


-1 


562 


19 ± 43 






E3.2 


Ks 


-3.093 


-74 ± 38 


-1 


525 


84 ± 24 






E3.3 


Ks 


-3.252 


-262 ± 1 / 


-1 


491 


2 it 15 






E3.4 


Ks 


-3.203 


-97 it 20 


-1 


632 


-23 i 27 






E3.5 


TT 

n 


-3.026 


-12 it 51 


-1 


501 


-126 it 48 






E4.0 


H 


-3.210 


-227 ± 4 


-1 


409 


26 ± 4 


80 ± 70 


-30 ± 80 


E4.1 


Ks 


-3.281 


-211 ± 21 


-1 


410 


206 ± 18 






E4.2 


H 


-3.364 


-25 it 43 


-1 


402 


-80 ± 27 






E5.0 


Ks 


-3.382 


-136 ± 12 


-1 


502 


175 ± 14 






E5.1 


Ks 


-3.391 


-87 it 38 


-1 


569 


42 ± 34 






E5.2 


Ks 


-3.474 


-340 ± 38 


-1 


541 


3 it 28 






E5.3 


Ks 


-3.458 


-139 ± 39 


-1 


607 


163 ± 30 






E5.4 


Ks 


-3.446 


-68 it 38 


-1 


469 


30 ± 34 






E6 


H and Ks 


-3.010 


-133 ± 7 


-1 


423 


23 ± 7 






E7 


Ks 


-3.495 


-227 ± 29 


-1 


391 


167 ± 28 






E8 


Ks 


-3.320 


-118 ± 18 


-1 


866 


95 ± 20 







Note. — The positions refer to reference epoch 12-05-2005 and arc measured relative to Sgr A*. 
The dominating uncertainty of the position is caused by the coordinate uncertainty of about 2 
mas jReid et al]l2007l : ITrippe et"alll2008l : IGillessen et al.ll2009l ') . 



would be detectable against the narrow emission hues. 
Given that a WR-star cannot be fainter than rnxs = 16.5 
(jCrowthcr 2007), only a few types of WR-stars could be 
hidden in E3.0. Also, E3.0 does not resemble a dusty 
WR-star (of which E4.0 is one), in which case the lines 
would be diluted but still visible. 

It is worth noting that also the background emission 
is well determined from our data. The seeing halo of 
E3.0 is of no concern since the lines in the background 
are stronger in relative terms. The extinction appears 
to be similar in the whole region, since we do not mea- 
sure any change in the relative line strengths for lines 
from one ionization state of an element over the spectral 
cube. In contrast, the line strength ratios of different el- 
ements differ between the background and E3.0 (see for 
example the HI, Fell pseudo doublet at 1.64 /im in Fig- 
ure 01). These spatial changes in the spectral properties 
appear to be gradual, consistent with a gaseous nature 
of E3.0. They are strongest at the position of E3.0, but 
they do not correlate in general with the positions of the 
red sources in IRS13E. 

The radial velocity of E3.0 determined from the 
Brackett-7 line is wlsr ~ — 25km/s, n ieasured at the 
ositio n of the continuum peak of E3.0. iPaumard et al.l 
20061) measured a velocity of wlsr = 87 ± 20km/s for 
E3 but did no t iden tify the type of star. Probably 
IPaumard et al.l ()2006D measured an average of E2 and 
E4.0 instead of the E3 velocity. 

Given the photometric variability of E3.0, we fit the 
SED based on the narrow-band data set from 2004 and 
the one based on the SINFONI spectroscopy from 2009 
separately. Using the latter we obtain Aks = 3.89 ±0.16 
and a warm dust component of T = 970 ± 17 K, see 
table [2] and figure [51 The la flux limit on the stellar 
component is mxs = 16. Using the earlier imaging data 
yields consistent results. Also in this case, no stellar 
component is needed, however the limit is less constrain- 
ing due to the smaller wavelength coverage. The extinc- 
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Fig. 4. — H-tK-band spectrum (red) of the brightest component 
of IRS13E3: E3.0. For comparison the spectrum of the nebular 
background is shown, too (blue, shifted). In the upper part, the 
line positions for various transitions are marked. The weak CO 
band heads in the background spectrum are caused by stray light 
of late-type stars. The weak absorption feature of E3.0 around 
Brackett 7 and Helium are likely caused by the slightly imperfect 
subtraction of the stellar emission lines. 



tion value = 3.5 ± 0.7 is consistent with the former 
one. We also obtain an independent extinction estimate 
of Aks = 3.44 ± 0.15 from the Bra ckett series of E3.0 
by using the relative line strengths of iHummer fc StorevI 
(^198'^) and the el ectron densities and temperatures from 
iLutz et al.l (|1993[ ). Our different extinction estimates ap- 
pear to be sufficiently consistent with each other. Using 
a value around ^ks = 3.4 in conjunction with the 2009 
data yields that the magnitude of a possible stellar com- 
ponent is less than toks = 18. 
We also compare the spectrum of E3.0 with those 
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Fig. 5.— SED data and model fits for El, E3.0, E3.5, E4.1, E5.0, 
E8, the minispiral and the sum of E3.1 and E3.2. All data and fits 
are offset from each other for clarity. For all sources the extinction 
correction of El is applied. El and E3.5 are fit well by single 
blackbodies of stellar temperature. The objects E3.0, E3.2+E3.3, 
E4.1, E5.0 and E8 can be fit as single blackbodies with a dust- 
like temperature. Stellar components do not improve these fits. 
The SED of the minispiral shows two blackbodies: A warm dust 
component and stellar stray light. 



of young stellar objects (YSO). Such objects are not 
expected to reside in the central parsec, because the 
gas density t here is much t o o low for star forma - 
tion ( C hristo pher et all [20051: iBonnell fc Ricd I2008D . 
lEckart et al.i (j2004D and IMuzIc et alJ (|2008l ) speculate 
that some of the objects of IRS13N nevertheless could 
be YSO, arguing that shocks could allow star formation 
to take place. The combination of color and luminosity 
for E3.0 argues against a typical Herbig-Be-star. Also, 
no distinctive feature of a massive YSO is visible in E3.0. 
Except for the Brackett lines, the lines observed in E3.0 
are unknown (e.g. Felll) or weak (e.g. He) in massive 
YSO. A massive YSO does not necessarily show other 
lines apart from Brackett-7 in the K-band, but often CO- 
band h eads, H2, or the Pfund series in emission can b e 
visible ()Bik et all [20061 iMartin- Hernandez et al.ll2008D . 
Such a clear indication for a YSO is not observed in 
E3.0, nor in other parts of IRS13E. Furthermore, there is 
no additional local extinction towards E3.0 while YSOs 
show locally enhanced extinction. We conclude that it is 
unlikely that E3.0 is a YSO. 

The source E3.0 in the H-band has a velocity simi- 
lar to the the radio source in IRS13E (,Zhao et al.ii2009i ). 



Therefore the H-band emission is likely due to dust as- 
sociated with the gas, which causes the radio continuum 
emission. Hence, also qualitatively a higher extinction 
value (like for example expected in YSO) seems not to 

be needed. 

iMaillard et al.l (|2004[ ) found two similarly bright com- 
ponents E3A and E3B, which is different from our result 
at higher resolution. Their additional stellar components 
do not fit our SED for E3.0. Our flux at 1.48 ^m is a 
factor « 1.8 lower than the sum of the flux models for 
E3A and E3B, but at 2.18 //m our flux is only a fac- 
tor R:^ 1.15 low e r. Th e hint for the stellar component in 
iMaiflard et~all (|2004[ ) is only based on the 1.6 /im-flux 
measured in the filter F160W from NICMOS. Since this 
is a very broad filter, and E3.0 is a very red source, the 
central wavelength of the filter is considerably bluer than 
the effective wavelength. Hence, using the central wave- 
length will overestimate the fiux. We speculate that this 
effect made the stellar component necessary. Also, the 
amount of stellar flux contributing at that wavelength to 
the total flux is small, and hence a smal l measurement 
error might change the conclusion from I Maillard et alJ 

iooi. 

To summarize; the potential stellar component in E3.0 
is fainter than most known WR-stars. A main sequence 
star in E3.0 would need to be fainter than B1.5V. This 
seems unlikely, however, since these stars do not produce 
dust. With three different methods we obtained consis- 
tent extinction values, similar to the extinction towards 
El and E2. We conclude that E3.0 is neither a star nor a 
YSO, but rather a concentration of warm dust and gas. 

The source E3.0 is not the only dust source of its 
brightness in the GC. For example, the source IRS2L re- 
sembles E3.0 in many respects and also shows extended 
H-band emission. 

4.2.2. E7, E8 and the surrounding gas 

The background around IRS13E is not only visible in 
gaseous emission lines (figure |4]) but also in the contin- 
uum. The deconvolution artifacts hamper any photomet- 
ric measurement of this continuum. Instead, we use the 
objects E7 and E8 (figure[2|) that are less affected. These 
objects are farther away from the center of IRS13E than 
the other sources and likely do not belong to IRS13E but 
to the background. Spectroscopically, E8 appears to be 
indeed a local overdensity in the gas. 

E7 and E8 are not detected in the H-band. We fit a 
blackbody model to each of the two objects, fixing the 
extinction to its upper bound at ^ks = 3.9. We find 
that no star is needed to explain the SED of these two 
objects. We obtain a temperature of T = 692 ± 55 K for 
E7 and of T = 732 ± 10 K for E8. Thus, E7 and E8 are 
colder than most objects in IRS13E. 

We also detect three faint stars (mxs ~ 15) in the 
SINFONI data at the approximately same distance (w 
450 mas) to E3 .0 . All these stars show CO band heads. 
iBuchholz et al.l (|2009f ) correctly identify the star with 
the strongest CO-band heads as late-type, the others 
are misidentified as early-type stars. The radial veloc- 
ities of these stars are different from the radial velocity 
of IRS13E. Thus, the stars differ from IRS13E not only 
in proper motion (jTrippe et al.l 120081) but also in stellar 
type and radial velocity. We conclude that they are fore- 
or background objects. 
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The bac kground around IRS13E is part of the 
minispiral (iLo fc ClaussenI 1 19831 iPaumard et alj 120041: 
I Zhao et "all I2OO90 . The minispiral is visible in our L'- 
and Ks-band data. In order to measure its SED, we se- 
lect seven boxes of w 0.2" diameter far away from bright 
stars at a distance of around 1.5" to IRS13E. We exclude 
the narrow-band data around Brackett-7 because of the 
strong emission line there. We calibrate the flux from 
the narrow-band filters by fitting these flux values as a 
linear function of wavelength and setting equal the flux 
at 2.18 /im to the measured Ks-band flux. We fit two 
blackbodies to the SED data (figure E]) . One accounts 
for the stellar stray light, the other represents the dust. 
The fits for all fields are similar. The average extinction 
is Aks = 3.64±0.55±0.56, consistent with the extinction 
value derived for El. The average dust temperature is 
T = 585 ± 22 ± 33 K. Thus, the minispiral is colder than 
the dust in IRS13E. 

One of th e seven fields conta i ns the central sources 
of IRS13N (lEckart et al.l [2Q0l IMoultaka etaH 120051: 
iMuzic et al.l I2008D . In this case the fit yields ^ks = 
3.7 ± 1.6 and T = 664 ± 174 K. This argues against a 
highly extincted object (like a YSO) as these authors 
had proposed. The object IRS13N/3 is within the field of 
view of our SINFONI data. We do not find any signifi- 
cant differences between its spectrum and the spectrum 
of the gas around IRS13E. Also, the velocity of IRS13N 
in L'-b and (Muzic ct al. 2008!) and at radio wavelengths 
(jZhao et a l. 2009) make it more likely that IRS13N ac- 
tually consists of dust blobs. They appear denser than 
in other regions in the GC, but still they are much less 
prominent th an E3.0. The com mon proper motion of 
these sources (|Muzic et al.ll2008() could also be explained 
by a recent formation of the dust clouds. 

5. PROBABILITY FOR AN IMBH IN IRS13E 

We now discuss the nature of IRS13E considering our 
NIR observations. There are four main explanations: 

• IRS13E is a cluster in the process of dissolu- 
tion, one of variou s possibilities suggested by 
ISchodel et aLl (pfiol . This is very unlikely, be- 
cause it would get diluted extremely quickly. In the 
case of IRS1 3 N wit h a similar dispersion as IRS13E 
IMuzic et al.l ()2008D show in their figure 2 that such 
an association is diluted quickly in the GC. We 
estimate that after just 200 years one would not 
recognize anymore the overdensity of three young, 
massive stars of IRS13E as such. We exclude this 
option. 

• IRS13E is a cluster bound by stellar mass. From 
our detection limits, we estimate that the total stel- 
lar mass of the cluster could at most be 2000 Mq. 

• IRS13E is a cluster bound by additional mass, e.g. 
an IMBH. 

• IRS13E is a chance association. 

5.1. Stellar surface density in IRS13E 

A physical cluster leads to an overdensity in projection 
compared to the local background density. Detecting a 
local, significant overdensity in the stellar surface den- 
sity would hence be a strong indication that IRS13E is a 
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Fig. 6. — Surface densities of objects in the core (r = 0.27") of 
IRS13E. The errors are Poisson errors. The upper limits are Icr- 
limits. Red: All objects. Blue: Cluster member candidates only. 
The red and blue data are offset slightly in the horizontal direction 
for clarity. The solid line shows the estimated background density 
for our data. 

cluster. From our data we meas ure the projected space 
densities of stars in IRS13E. Pau mard et al.l (120061 ) found 
13 objects with tor < 19.4 inside a radius of 0.3" in 
IRS13E. Similarly, we find 15 objects in H- and Ks-band 
with mn < 18.5 inside a radius of 0.27". 

We determine the density of background stars in a non- 
confused region with a size of 5.15 sqas north of IRS13E 
for stars fainter than mn = 14. Because there the ex- 
tinction is higher than towards I RS 13 E according to the 
extinction map of lSchodel et"all ()2009t ) , we shift the mag- 
nitudes of all the stars in this field by 0.4 magnitudes to 
brighter magnitudes. F or stars brighter th an mn = 14 
we use the results from iGenzel et all ()2003D to estimate 
the background density. We find that the density of ob- 
jects in IRS13E is higher than in the surrounding field 
at nearly all magnitudes (figure [6]), if all objects are con- 
sidered. However, only for six of these objects a star 
is required to explain the observed SED. For the other 
objects, a possible stellar component needs to be fainter 
than about ttiks ~ 15.5. This corresponds to a stellar 
surface density of 26 as~^, which is five times higher than 
the background. However, E6 (figure [Tj) is not a member 
of the cluster, given its age well above 6 Myr. Further- 
more, the two faint stars E3.5 and E4.2 have magnitudes 
consistent with the red clump and hence are likely also 
older than the bright stars. In the magnitude bin of E3.5 
and E4.2, the background population accounts already 
for one of the two stars. Thus, the overdensity in this 
bin is insignificant, too. 

We conclude that the magnitude integrated overden- 
sity of objects in IRS13E (with a chance probability of 
2 X 1 ~ ^ , IPaumard et all (|200(l ) ) is caused by dust clumps 
and not by stars. This overdensity cannot be used as ev- 
idence for a cluster. Nevertheless, the concentration of 
the three bright early-type stars in IRS13E remains puz- 
zling, an argument which we analyze quantitatively in 
section I5T3I 

5.2. The structure of IRS13E 

A significant number of dust clumps in t he GC are 
found far away from stars (|Muzic et al.ll2007j) . But most 
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TABLE 2 

Spectral properties of the objects in IRS13E 



Name 


data 


spectral type 


Aks 


Ti \K] 


To fKl 

L J 


star mj ks 


dust 7722 Ks 


El 


NB, H, K, L 


OB I 


3.64 ± 0.07 


26000 




10.38 ± 0.06 




E2 


SINF, NB, H, K, L 


WN 8 


3.51 


30000 


712 


10.43 


14.74 


E3.0 


SINF, NB, H, K, L 


gas 


3.89 ± 0.16 


9480 


970 ± 17 


> 16 


11.42 ± 0.05 


E3.1+E3.2 


NB, H, K, L 


gas 


3.9 


9480 


955 ± 24 


> 18.03 


13.24 ± 0.18 


E3.3 


NB, H, K, L 


gas 


3.9 


9480 


987 ± 54 


> 15.16 


13.65 ± 0.35 


E3.4 


NB, H, K, L 


gas 


3.9 


9480 


954 ± 49 


> 16.53 


14.15 ± 0.36 


E3.5 


NB, H, K 


3.49 ± 0.36 


4000-19000 




15.30 ± 0.26 




E4.0 


SINF, NB, H, K, L 


WC 9 


3.64 


41000 


1400 


11.95 


12.51 


E4.1 


NB, H, K, L 


gas 


3.9 


9480 


844 ± 24 


> 18.36 


14.38 ± 0.21 


E4.2 


NB, K, H 




3.9 


4000-19000 


517 ± 252 


15.77 ± 0.11 


16.64 ± 2.5 


E4.3 


K 












15.48 ± 0.10 


E5.0 


SINF, NB, H, K, L 


gas 


3.4 


9480 


700 ± 8 


> 18.46 


14.08 ± 0.10 


E5.1 


NB, H, K, L 


gas 


3.9 


9480 


947 ± 25 


> 17.75 


14.25 ± 0.14 


E5.2 


NB, K, L 


gas 


3.9 


9480 


642 ± 70 


> 15.6 


15.49 ± 0.10 


E5.3 


NB, H, K 


gas 


3.9 


9480 


978 ± 524 


> 16.04 


15.1 ± 1.6 


E5.4 


NB, K 


gas 


3.4 




817 ± 474 




15.54 ± 2 


E6 


NB, H 


K3III 


3.68 ± 0.09 


4300 




13.82 ± 0.10 




E7 


NB, K, L 




3.9 


9480 


692 ± 55 


> 16.71 


16.02 ± 0.40 


E8 


NB, K, L 


gas 


3.9 


9480 


732 ± 10 


> 17.14 


14.51 ± 0.11 



Note. — We fit one or two blackbodies with temperatures Ti and T2 and one extinction parameter to our SED data. We 
allow extinction values between Aks — 3.4 and Aks — 3.9 and clip it if needed at the respective edge of the interval. We derive 
the temperature of the stellar component from an atmosphere model or from the observed lines when possible. Otherwise the 
temperature is set to T — 9480 K. All limits are iCT-limits. In column two, SINF means the SED fitting results are based on 
SINFONI data, NB that the narrow band images arc used. 



of these dust clump are invisible in Ks-band and none of 
them is visible in H-band. In the H-band the dust blobs 
in IRS13E are among the brightest non-stellar structures 
in the central parsec. The dust clumps close to E3.0 
are nearly twice as hot as the minispiral and also denser 
(« 4x for E3.0, estimated from surface brightness and 
temperature) than the minispiral. The brightest and 
hottest dust blobs (around E3.0) are located between 
the two bright WR-stars E2 and E4.0. The star E2 is 
the brightest of all WR stars in the GC and E4.0 is a 
dusty WR with the hottest dust in IRS13E. 

Therefore E4.0 might be the most important dust 
source in IRS13E. But it seems unlikely that E4.0 is re- 
sponsible alone for the bright dust blob E3.0, given that 
it does not move away from E4.0, while the motions of 
the blobs E4.1 and E4.3 are consistent with being out- 
flows from E4.0. Instead, we assume that the winds of 
E2 and E4.0 collide in the E3 region and cause the bright 
dust blobs in this way. For this model it is not necessary 
that the dust originates from the stars but only that it 
is heated by the two stars. The fact that the velocities 
of these blobs are not pointing towards either of the two 
WR stars supports this picture. Also, the eastward com- 
ponent of the R.A. motion of the blobs increases from 
west to east as expected in this model. Accordingly, E2 
and E4.0 would not only be close in projection, but would 
indeed have a small 3D distance. 

A potential weakness of that model is that the motions 
of the blobs E3.3 and E4.1 do not seem to fit. But this 
is only a worry if no other forces than repulsion from 
the stars act on the winds. For example, the presence of 
the gas from the minispiral might decelerate the winds. 
A second problem might be that it is not clear whether 
the conditions for dust formation are satisfied, give n that 
our s ystem is very different from WR140 (Monnic r et al.l 
[2OOI . 

Most likely El does not contribute to the dust produc- 
tion. It is further away and supergiants have weaker 
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Fig. 7. — Velocities of the objects in IRS13E relative to the mean 
motion of the two WR-stars E2 and E4.0. The near-infrared mo- 
tions are marked with black arrows, the boxes denote the 1-cr er- 
rors. Stars are marked with asterisks (blue: young stars; green: 
stars without spectral identification). Dust clumps are marked by 
red disks. The sizes of symbol s indicate the Ks- band flux. The 
13 mm radio motions (violet) of I Zhao et al.l II2009I) are shown, too. 



winds than WR-stars. Furthermore, the velocities of 
some dust clumps are pointing towards El. Therefore 
El can, but needs not to be close in real space to the 
WR-stars. 

The dust clumps in E5, E7 and E8 fit less well into our 
simple model of two WR stars because of their velocities. 
They are also mostly colder than the dust clumps in E3 
and E4, but still warmer than the minispiral. In addi- 
tion, they are farther away from E2 and E4.0. Likely the 
minispiral is more important for their formation. 
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5.3. The overdensity of young, massive, early-type stars 

Because most of the fainter objects in IRS13E are dust 
clumps, the most prominent feature of the potential clus- 
ter IRS13E is the concentration of the three bright early- 
type stars. Therefore we quantit atively evaluate the ar- 
gument of iMaillard et al.l (|2004[ ). that this association 
cannot be a chance occurrence. On the one side, the 
number of stars is very small, and one could argue that 
IRS13E cannot be a significant overdensity. On the other 
hand, five or even six phase space coordinates of these 
stars are similar. 

We calculate the probability that a star association like 
IRS13E occurs randomly, i.e. a pseudocluster, by means 
of Monte Carlo simulations. IRS 13E cannot be part o f 
the clockwise disk of young stars (jPaumard et al.ll2006l ). 
Consequently, the size of the sample of stars considered is 
chosen to equal the observed population of those O- and 
WR-stars that are not part of the clockwise disk, which 
yields 96 stars with r < 10". We assume that as many of 
the counterclockwise moving sta rs are in a disk-like con- 
figuration as suggested by Bartko et al.l (2009) ; the rest 
is on randomly oriented orbits. The disk thickness and 
eccentricity distribution of the counterclockwise disk are 
not well known. The thickness is likely larger than 10° 
(the value for the clockwise disk) and smaller th an 25°, 
becau se a thicker disk would be invisible (Bartk o~et al.l 
I2009D . The eccentricity distribution is more uncertain. 
For each realization we have a set of ecc entricities with 
an rms width of 0.22 (jBartko et al.|[2009l) and a median 
eccentricity which is varied for each realization. The me- 
dian values are chosen such that they are distributed 
symmetrically around the mean eccentric ity of the clock- 
wise disk, in the range fr om 0.17 to 0.60 (jPaumard et al.l 
l2q06HBartko et al.l[2009l ). We use these eccentricity dis- 
tributions for the counterclockwise disk and for the ran- 
dom stars. The impact of the width of the eccentricity 
distribution on the derived probabilities is small com- 
pared to influence of the median eccentricity. The distri- 
bution of semi-major axes is chosen such that the final 
surface density profile resembles the one observed for the 
stars that are not part of the clockwise disk. A change 
of the semi-major axis distribution compatible with the 
errors of the observations has nearly no infiuence on the 
derived probabilities. The actual number of stars is var- 
ied, too. On purpose, we don't restrict the sample to the 
stars at least as bright as the IRS13E members. Fainter 
0-stars are barely less massive, and dynamically they 
would not behave any differently (jPortegies Zwart et al.l 
[200l . 

The criteria for finding a cluster from the simu- 
lated stars are motivated by the observed properties of 
IRS13E. If anywhere in the simulated set three stars are 
found for which the following four conditions are met, 
the set is marked as containing a pseudo-cluster of three 
stars. The following cut values are the measured values 
plus the 1(7 errors: We require that rdustor, 2D ^ 0.19" 
and that the proper motion dispersion ctv, 2D < 94km/s. 
Furthermore, we demand that for any two stars of the 
three r^u < 0.46" and for any (other) two stars of them 
Avz < 76 km/s, A pseudo-cluster of two stars is de- 
fined by the criteria d2D < 0.33", i7v,2D £ 23 km/s, and 
dsD < 0.46". In the following we give the probabilities 
to find such configurations, and in brackets we give the 
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Fig. 8. — Probability for observing a cluster of two (red) or three 
(blur) early-type stars that resemble IRS13E in a background pop- 
ulation of given sample size. The errors are derived from assuming 
different distributions of orbital elements for the simulated stars. 
The red and blue data are shifted for clarity horizontally. The back- 
ground population consists of the O/WR early- type stars outside 
of the clockwise disk. The violet bar marks the range applicable 
for the GC. 



values if the requirement that two stars need to be close 
in the ^-coordinate is omitted. 

By randomly drawing a stellar distribution many times 
and varying the assumptions on the dynamical properties 
of the O/WR stars population, we find that a triple of 
stars matching the criteria for a pseudo-cluster of three 
stars out of 96 has a chance probability of 0.2 % to 2% 
(0.6% to 4%) with the most likely probability being 0.8%. 
Finding a pair of stars out of 96 that matches the crite- 
ria for a pseudo-cluster has a chance probability of 0.4% 
to 5% (4% to 14%, figure El). Demanding in addition 
that E2 and E4.0 have similar masses lowers the chance 
probability for such a pair to 0.1%, since only 28 of the 
96 stars from which we draw the pseudo-clusters have 
masses comparable to that of E2 and E4.0. 

The dominant error source here are the assumptions on 
the dynamical properties of the O / WR stars population, 
as indicated by the ranges of values in the previous para- 
graph. The uncertainties of the probabilities induced by 
the measurement errors of the proper motions are very 
small. They lead to relative errors of the probabilities 
well below 10%. Also the influence of Sv^ on the prob- 
ability is small with 13% (relative error). We conclude 
that the occurrence of a pseudo-cluster similar to IRS13E 
is not expected, but still not impossible^. 

We also note that E2 and E4.0 could be a binary in 
an orbit with a pericenter distance of more than 1.3 pc 
(see section l5.4.2|) . Heavy bi naries indeed often consist o f 
two similarly massive stars (jKobulnickv fc Frved[2()?)l . 
On the other hand, the hypothesis that all three stars 
are physically connected and reside outside the central 
parsec seems improbable since such young, massive stars 
actually are very rare in the Galaxy. And regardless of 
the location of IRS13E along the line-of-sight, the high 
velocity dispersion would require additional, unseen mass 
to bind the cluster. 

* In the clockwise di sk one also find s a comoving group of stars: 
The IRS16-SW group ijLu et al.ll2005I V This group is less concen- 
trated than IRS13E. 
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5.4. Cluster mass from velocity dispersion 

Because the probability for IRS13E being a chance as- 
sociation is relatively small, it is reasonable to exploit the 
idea that IRS13E indeed is a cluster. Therefore we esti- 
mate the mass needed to bind the cluster by analyzing 
the velocity dispersion. 

5.4.1. Assuming El, E2 & E4-0 are cluster members 

Classical mass estimators are very uncertain when ap- 
plied to a sample of three stars only. Hence, we simulate 
clusters of three stars around a central mass in a Monte 
Carlo fashion for deriving the probability distribution for 
the value of the mass. For this one needs to marginalize 
over many parameters, mainly the 18 orbital elements of 
the three stars. Since we assume here that IRS13E is a 
cluster, we normalize the probability such that the total 
probability of clusters accepted over the full mass range 
is 1. 

We assume that all angles and phases are randomly 
distributed and that the eccentricities follow a thermal 
distribution n(e) oc e. The semi-major axes a are drawn 
from a = 0.4" xp*^'^, where p is a random number in [0,1]. 
This makes the distribution of projected distances similar 
to the observed one. In this way we simulate 5 x 10* 
clusters around different IMBHs with masses between 
300 M0 and 3 x 10^ Mq in steps of AM/M = 0.15. 

The peculiar property of IRS13E is that one star has a 
significantly different velocity compared to the other two. 
Thus, we choose from the simulated clusters the ones, 
which have nearly the same velocity differences as the 
observed cluster: Between any two stars this is | V2D. sim — 
^2D,obs| < 10 km/s and between two of them additionally 
|V^,sini — Ki^obsl < 39 km/s. We neglect the exact spatial 
distribution of the stars in IRS13E. But on average the 
simulated stars from the chosen clusters are at similar 
positions as the observed stars. 

The uncertainties of the derived probability density 
due to the measurement errors and due to the assump- 
tions on the semi-major axis distribution are small. 
Even though the exact shape of the distribution changes 
slightly, it will continue to be a concave function normal- 
ized to 1, with a broad maximum around 2 x 10* Mq. 

We derive that the most likely mass is Mimbh = 
16900 M0 with a 68.3 % confidence region from 4000 
to 59400 Mq. The probability distribution for A/imbh is 
highly non-Gaussian, see black curve in Figure [lUl The 
large uncertainty on the mass is due to the low number 
of stars. Nevertheless, the mass needed to bind IRS13E 
exceeds the mass of the three stars. Hence, if IRS13E is a 
cluster of three stars, it needs to be bound by additional 
mass that is not detected by our observations. This mass 
could be present in the form of an IMBH. 

5.4.2. Assuming only E2 & E4-0 are cluster members 

If one assumes that E2 and E4.0 form an equal-mass 
binary, the measured velocity difference yields a mass 
for each star of 180 Mq. The two stars have almost cer- 
tainly a much lower mass, but the argument illustrates 
that an IMBH with a mass similar to the equal-mass 
binary is sufficient to bind the two stars, too. Since 
this mass is less than the limit on the stellar mass in 
IRS13E (2000 Mq), one would not need an IMBH to ex- 
plain IRS13E. 



The simple binary hypothesis is problematic for other 
reasons. The tidal forces of the SMBH disrupt a 
binary with a separation larger than the Hill-radius 
ijBinnev fc TremaiM I2008D during pericenter passage. 
For a bound, equal mass binary there exists no value of 
the unknown line-of-sight distance \z\ that would make 
the binary survive a pericenter passage (figure [TT|). Also, 
the destruction timescale for such a soft binary in a dis- 
tance of 3.5" to the SMBH is only 19000 years (calcu- 
lated for th e observed projected separation of E2 and 
E4.0. iBiime v & Tremainc (2008)). Furthermore, a mass 
of 180 Mq per star seems unrealistic, but this mass esti- 
mate could be too high due to the measurement errors. 
A mass per star of 85 Mq is excluded only at the 2a- 
level. Taken together, these arguments make the binary 
hypothesis unlikely, and it seems more natural to assume 
that additional mass is needed for binding E2 and E4.0. 

5.5. Acceleration limits make an IMBH less likely 

A binding IMBH can manifest itself in the accelerations 
it imposes on surrounding stars. Vice versa, the absence 
of accelerations makes it more likely that IRS13E is a 
chance occurrence. Thus, we calculate now the proba- 
bility for observing the observed (insignificant) accelera- 
tions, firstly for the case that a binding IMBH is present 
and secondly for the case that no additional mass is 
present in IRS13E. 

5.5.1. Binding IMBH present 

From our data we can only test for accelerations in the 
proper motions: 

r2D G AflMBH ^ , , r2T) 
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R. 



3D 
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IMBH 



(in + ^^r^' 
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We do not measure any significant accelerations (ta- 
ble [T|). Firstly, this sets only limits on the location of 
the stars relative to the IMBH along the line of sight (z). 
But secondly, it is unlikely that a cluster extends much 
more in the z-dimension than in the plane of sky, and 
this actually constrains the mass A/jmbh of the IMBH. 

We use again as in section section [5^ the Monte Carlo 
simulated cluster, which are similar to IRS13E. For all 
this clusters, we calculate the expected values of 020- 

Figure ini shows that E2 is nearly not constraining, since 
for nearly all simulated clusters the expected value of a2D 
is smaller than the la acceleration limit. In case of E4.0 
the expected accelerations are similar to E2. Since the 
measurement errors of E4.0 are four times larger than 
for El and E2, this star is of least importance. As con- 
sequence, only the acceleration limit of El is important. 
This is due to the fact that its velocity differs fro m the 
one o f E2 and E4.0. As pointed out by Schodcl et al.l 
([2005) in their figure 2, the IMBH is likely closest to El. 
Therefore its motion is most sensitive to the mass of the 
IMBH. 

It is apparent in figure [9] that El excludes an IMBH 
of lower mass more efficiently than an IMBH of higher 
mass. This seems paradoxical, but results from the fact 
that the lower the mass is, the more likely it is that the 
IMBH resides close to El. If the IMBH is close to El, the 
full velocity difference between El and the other stars is 
due to El and one can write. 
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Fig. 9. — Plane of sky accelerations of the simulated stars in 
our Monte Carlo clusters. Every red dot marks an acceleration of 
El in one of the 50000 simulated clusters. The blue dots stands 
for E2 in the same simulations. The horizontal line represents our 
acceleration limit. The region above the line is excluded at the 
la-level. 



Reducing Mimbh thus lowers tei if the velocity difference 
has to stay at the measured value. For the acceleration 
oei follows 



A^IMBH 
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This anti-proportionality of acceleration to Mimbh is ap- 
parent in the simulations (red points in figure [9]). 

We now calculate how probable the observed accelera- 
tions are, if one assumes a certain M imbh ■ Hence, we use 
those simulated clusters from section lS^ that fall into the 
corresponding mass bin and calculate for each case the 
for the observed accelerations given their errors and 
the simulated acceleration values. This in turn is then 
converted to a probability. Because all simulated clusters 
are equally likely, we use the average over this mass bin 
as probability for obtaining the measured accelerations 

(figure nni). 

All IMBHs with Mimbh < 10000 M© have a probabil- 
ity of less than 0.9% to be consistent with our accelera- 
tion limits. At the most probable mass of 14000 Mq the 
probability of not detecting an acceleration is 2.5%. For 
all mass bins the probability of Mimbh > is less than 
3%. The mass that yields acceleration values that are 
least unlikely is about 60000 Mq. 

The main uncertainty for the derived probabilities is 
due to the assumption on the semi-major axis distribu- 
tion a = 0.4" X Using 0.28" or 0.52" instead of 0.4" 
(which makes the simulated clusters still look reasonably 
similar to IRS13E) changes the resulting probabilities by 
up to a factor of « 2. 

5.5.2. No additional mass 

If there is no IMBH in IR13E we expect zero acceler- 
ation. We calculate the to the accelerations for this 
case, too, and obtain a probability of 39 % for observing 
the measured acceleration values under the assumption 
that no additional mass is present in IRS13E. 

The mass estimate of w 1000 Mq from the assumption 
that only E2 and E4 are bound would yield an accel- 



eration of w 3/ias/yr^. This is much smaller than our 
measurement accuracy and we cannot test this hypothe- 
sis. 

5.6. Proper motion of Sgr A * 
5.6.1. Binding IMBH present 

iGualandris fc Merrittl (|2009f ) show that the proper mo- 
tion of Sgr A* perpendicu lar to the plane of the Galaxy 
(|Reid fc Brunthaleii[200l is the most promising method 
for excluding an IMBH at distances similar to that of 
IRS13E. For IRS13E itself, we can take advantage of the 
fact that we know the velocity of the potential cluster 
and hence of the proposed IMBH. Thus we can predict 
the expected reflex motion of Sgr A* via momentum con- 
servation. We do n ot need to assume circular m otion in 
the plane of sky, as IGual andris fc Merritti ('2009') had to. 

We use the value from Rcid fc Brunthal cr (2004) who 
find that the motion perpendicular to the Galactic plane 
is —0.4 ± 0.9km/s. We estimate the uncertainty of the 
motion of the IMBH simply by the velocity dispersion of 
the potential cluster IRS13E. 

We proceed as in section 15.51 For a given mass bin, 
we calculate the expected reflex motion of Sgr A* due 
to the IMBH the motion of which is varied by the esti- 
mated uncertainty. We compare the such calculated re- 
flex motion with the measured value. The corresponding 
probabilities we again obtain from the values. Finally 
we use the average per mass bin to estimate the prob- 
ability of the radio data of Sgr A* under the assump- 
tion that an IMBH of that mass is present (figure [T0| . 
The small insignificant proper motion of Sgr A* is quite 
probable for small (< 2OOOOM0) IMBH masses. But for 
Mimbh > 46OOOM0 the data are inconsistent with the 
model at the 90% level. 

5.6.2. No additional mass 

We also calculate probability for that the measured 
motion occurs when no additional mass is present. We 
obtain a probability of 65.7%. 

5.7. Total probabilities 
5.7.1. Binding IMBH present 

In section [5.4.1l we have shown that if an IMBH massive 
enough to bind IRS13E is present, its most likely mass is 
(given the velocity differences between El, E2 and E4) 
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The corresponding distribution of probabilities 
Ppos,voi(w) is the solid, black curve in figure fTOl 
which is normalized to 1, since it assumes the presence 
of an IMBH. 

Next, we calculate combined probabilities. For each 
mass bin m, the probability that the measured acceler- 
ation limits and the radio data are compatible with an 
IMBH of mass m is 



Pace, radio 

(m) 

— Pace 

(m) X _Pradio(TO) 



(5) 



(figure [TOl red, dashed curve). Hence, at any given mass 
the probability for the data under the assumption that 
an IMBH is present is smaller than 0.8%. 

Further, we can calculate the most likely probability 
value, by evaluating a weighted average of Pacc, radio (w) 
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Fig. 10.— Probabilities for different IMBH masses, if IRS13E 
is a gravitationally bound star cluster of 3 stars (E1,E2, E4.0). 
Solid lines are probability densities per logarithmic mass bin of 
size AM/M = 0.15 and are normalized to 1. Dashed curves are 
probabilities calculated binwise. The thick regions mark the re- 
spective 68.3 % confidence regions of the most likely masses. Black, 
solid: probability distribution of the mass of the IMBH from the 
measured velocity differences. Blue, dashed: probability to ob- 
serve the measured acceleration limits for an assumed mass of the 
IMBH. Magenta, dashed: probability to observe measured radio 
motion of Sgr A* for an assumed mass of the IMBH. Red, dashed: 
combining acceleration and radio information yields for every mass 
bin the probability that the observed values for the accelerations 
and radio motions occur. Green, solid: Probability distribution of 
the mass of the IMBH, combining velocity, acceleration and radio- 
based information, 



over all mass bins, using as weights Ppos,vci(^), i-e. 

/ Pace, radio (jTl) X Ppos, vol(?w) dm 



Pcomb 



/ Ppos,voi("i) dm 



(6) 



where the denominator equals 1 due to our normaliza- 
tion. The result is 



Pcomb = 0.4%. 



(7) 



This value is uncertain by a factor of ~ 2. 
Normalizing the nominator from equation [5] to 1, i.e. 

/ Pace, radio (.Tlj 'Pnas^vcXXJ^) /„n 

Peomb(m) = — — ^- , , , (8) 

J Pace, radio ("ij X Ppos,vel(TO) dm 

is equivalent to a Bayesian approach, for obtaining a com- 
bined probability distribution of the mass. This results 
in the green, solid curve in figure 1101 It is normalized to 
1 by construction which reflects the fact that the pres- 
ence of an IMBH is assumed. The most likely mass of an 
IMBH is thus 



Ml 



comb 
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5.7.2. No additional mass 

For the case that no additional mass is present, we ob- 
tain Pace, radio (jTi = 0) = 26%. In scction 15.31 we showed 
that Ppos,vei('7i — 0) is 0.8%, which yields Peomb(m- = 
0) w 0.2%. The error on this value is a factor of w 3. 

6. DISCUSSION 

We are left with a puzzle. The probability to obtain 
the observed data (positions, velocities, accelerations, ra- 
dio reflex motion of Sgr A*) when no IMBH is present 
is 0.2%. Similarly unlikely is the probability of getting 




z [as] 

Fig. 11. — Periastron distance of IRS13E for possible locations 
along the line of sight (z) in the point mass potential of the SMBH. 
The dotted lines mark the Icr-error band. Considering tidal dis- 
ruption, only the z positions corresponding to periastron distances 
above the blue lines are possible for the respective masses. The two 
lower lines for 1000 and 20000 Mq assume a cluster with an IMBH, 
the upper line is for the case that E2 and E4.0 form a binary with 
150 Mq per star. The shaded region close to z = 0" corresponds to 
t he cases in which IR,S13E is in the counterclockwise disk (CCWS) 
(Bartko et al. 2009). The shaded region to the right corre sponds 
to a location of IRS13E in the minispiral j| Paumard et al.l 12004). 
which is possible for all configurations. It is not possible, that an 
IMBH binds IRS13E, and that at the same time IRS13E is located 
in counterclockwise disk. 



the same data if one assumes the presence of an IMBH. 
For any mass assumed, the probability is smaller than 
0.8%. Hence, both scenarios are similarly unlikely given 
their errors of about a factor 2.5. We therefore consider 
additional arguments. 

6.1. Dynamics of the young, massive stars 

The clockwise disk of early-type star s in the GC is 
significantly warped (jBartko et al.l [20091) . Such a warp 
can be the consequenc e of anoth er non-spherical mass 
((Navakshin 2005; Lock mann fc B aumgardt 2009), and 
thus an IMBH might manifest itself via the warp. How- 
ever, there are many other possibilities to explain the 
warp: The counterclockwise disk might suffice, or the 
warp can also happen in the gaseou s pha se during the 
disk formation (jHobbs & : Navakshin l2009() . Dynamical 
friction wil l make an IMBH spiral down to the SMBH 
(jFuiii et al., 2009). such that it would stay only for a very 
short time in the radial region of the disks. Hence, an 
IMBH is a less efficient warper per mass than a disk. 

One possible origin of a potential IMBH is the 
core c ollapse of a youn g, massive star cluster 
(|Porteg ics Zw art fc McMillan 2002). IRS13E could be 
the re mainder of such an ob ject that currently is spiral- 
mg m. ISchodel et all (I2005D firstly indicated th at usmg 
the simulations of IPortegies Zwart fc McMillanI ()2002[ ) a 
cluster of an unrealistic high mass is necessary for a such 
massive IMBH as in IRS13E. In addition, we compare the 
number of young, massive stars in the GC with the expec- 
tion for the most likely ma ss of a n IMBH of 20000 M^. 

In t he simulations of IPorteg ies Zwart & McMilla^ 
(|2002D : IPortegies Zwart et al.l ()2004. ): .Fuiii et al.. (,2000 
the IMBH obtains a mass of 0.1% to 9% of the total ini- 
tial mass of the cluster, depending on the assumptions 
on the IMF, the radii and winds of very massive stars, 
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the density of the cluster and whether a single super- 
nova stops the formation of the IMBH. Thus, for the 
most likely IMBH mass one needs an initial cluster mass 
of at least 2 x 10^ M© (jFuiii et alJ 12009). Assuming a 
Salpeter IMF between 1 and 140 M© (.Fuiii et alJl2009f ). 
stars more massive than 40 Mq would merge into the 
IMBH. The lighter stars should still be observable today, 
and one would expect to find 790 0-stars (20 — 40 Mq) 
and 1950 early B-stars (10 - 20 Af©). 

In the GC there are in total only w 290 0/WR-stars, 
correcting the currently detected number of such stars 
upward by a factor of 2.4, as estimated from incomplete 
areal coverage and detection incompleteness (figure 2 in 
iBartko et al.l ()2010D ). However, not all of them are com- 
patible with the idea that they originate from IRS13E. 

Figure [Tl] shows that it is not possible that an IMBH 
binds IRS13E, and that at the same time IRS13E is lo- 
cated in the counterclockwise disk. Hence, the inspi- 
raling cluster scenario requires that the orbits of the 
stars stripped off the cluster are found in edge-on or- 
bits in a common orbital plane. We estimate that at 
most 50 0/WR-stars could have such orbits. We ob- 
serve, however, a factor 15 fewer 0-stars than expected 
for an IMBH of 20000 M©. For the bright B-stars, the 
discrepancy is even larger. We conclude that the forma- 
tion of an IMBH of 20000 M© in IRS13E is inconsistent 
with the observed number of young stars. 

6.2. X-ray and radio emission 

Apart from Sgr A*, IRS13E is the only source with 
significant em ission detectable in the NIR, radio, and 
X-ray regime (|Zhao fc Gos£|[T998HWang et al.|[200l . At 
13 mm, the morphology of the minispiral is similar to the 
one in L'-band. This emission traces hot gas, probably 
associated with the dust visible in the N IR. 

The X-ray spectrum of IRS13E (Wan g et al.l [20061 ) is 
similar to the diffuse X-ray emission from th e local back- 
groun d and from the region around Sgr A*. iCoker et al.l 
(|200^ argue, that IRS13E is a wind nebula excited by a 
star, which left the Luminous Blue Variable (LBV) phase 
a short time ago. Maybe this scenario also works for the 
two WR-stars the winds of which probably collide in E3. 
The accuracy of the position of the X-ray source is not 
sufficient to decide whether the X-ray source can be iden- 
tified with any of the bright NIR sources in IRS13E. 

The observed flux in IRS13E is Lx(2 - lOkeV) = 
2 X lO^'^erg/s (jWang et al.l 120061 ) . corresponding to ~ 
10^^ Eddington luminosities (L^dd) for an IMBH of 
20000 M©. This number is much s maller than for th e 
X-ray binaries detected in the GC ()Muno et al.ll2005l ). 
The faintness of the X-ray sou rce in IRS13E was the 
main argument of Schodel et al.l (2005) against an IMBH 
in IRS13E. However, Sgr A* itself is radiating in qui- 
escence with only I 'x(2 - lOkeV) = 4 x 10~^^LEdd 
(jBaganoff et al.ll200lD such that a low radiative efficiency 
does not exclude the presence of an IMBH right away. A 
comparison on basis of the X-ray fiux only is difficult, 
since the environment and accretion state might differ 
significantly between Sgr A* and an IMBH in IRS13E. 

Th ere is no variability reported for IRS13E in the X- 
rays (|Muno et al.l 120051) , despite the fact that during 
roughly half of the observations of Sgr A* IRS13E was 
observed, too. Many black holes are variable in the X-ray 



domain (exa mples are X-ray binar ies (IMuno et al.ll2005D 
and Sgr A* ()Baganoff et al.ll2001l: iPorguet et a l. 2003)). 
In section [52] we argue that the dust seen in IRS13E is 
physically connected with the group of stars. Comparing 
the Br-7 emission in IRS13E and close to Sgr A* shows 
that the gas reservoir in the IRS13E is much richer than 
close to the SMBH. Taken together, this makes it unlikely 
that the reported, steady source is due to an IMBH. 

7. CONCLUSIONS AND SUMMARY 

We present a detailed analysis of the potential clus- 
ter IRS13E in the GC. We use AO-based images in H- 
and Ks-band from 2002 to 2009 for identifying objects 
in IRS13E and for performing astrometry. We detect 19 
objects on most high Strehl ratio images in Ks-band and 
15 of them also in H-band. The proper motions of the 
objects are well determined. 

In addition, we characterize the SED of the objects in 
IRS13E using AO-based integral field spectroscopy (1.45 
to 2.45 /im) for the brighter objects and (narrow-band) 
imaging from H- (1.66 /im) to L'-band (3.80 //m) for the 
fainter objects. We fit the SED of the objects by a model 
consisting of two blackbodies and an extinction param- 
eter. From this we conclude that 13 objects are dust 
clumps without embedded stars and that the extinction 
appears to be constant in the IRS13E region. Three ob- 
jects are fainter stars: One is a spectroscopically identi- 
fied late-type star. The two fainter stars of the magni- 
tude of the red clump are consistent with the expected 
number of stars per area at this magnitude bin. Thus, 
these three stars are likely background stars. Therefore, 
IRS13E is only a concentration of the remaining three 
early- type stars El, E2 and E4.0. 

They have the following properties: 

• E2 and E4.0 are probable located close in 3D, be- 
cause of the bright hot dust clumps between them. 

• A chance association of three early-type stars has 
a probability of around 0.2%. Depending on the 
assumptions on dynamical properties of the stel- 
lar population this value varies between 0.06% and 
0.6%. 

• In the case of a binding mass, we use Monte Carlo 
Simulations of a cluster of three stars around all 
reasonable masses of IMBHs (300 to 3 x lO^M©). 
From those simulations we select the ones that 
agree with the constraints from the observed po- 
sitions and velocities of the stars. The most 
likely mass is about 20000 M©. This mass exceeds 
the stellar mass (inclusive faint, invisible stars) in 
IRS13E. 

• We find that at any given mass less than 0.8% of 
the previously selected simulations agree with the 
observed acceleration limits of the three stars and 
the fact that Sgr A* appears to be at rest. This 
value is uncertain by a factor of w 2. 

• There are roughly 15 times too few young stars in 
the GC for the formation of an IMBH as massive 
as 20000 Mq. 

• The weak and non- variable X-ray source in IRS13E 
makes the presence of an IMBH unlikely. 
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Overall, we conclude that it is more likely that IRS13E does not host an IMBH. 
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APPENDIX 

A: EXTRACTING THE SPECTRAL ENERGY DISTRIBUTIONS 
Deconvolution of SINFONI data 

In order to distinguish between objects in the SINFONI data, we deconvolve the cube in smaU bandpasses, taking 
into account the variation of the PSF with wavelength. We extract line maps at three wavelengths at which single 
stars prominently dominate the respective images: at 1.49 //m the dominant star is E4, at 1.70 /im and 2.11 /xm it is 
E2. We manually remove apparent artifacts and smooth the PSF wings azimuthally. We then spectrally inter- and 
extrapolate these single- wavelength PSFs over the whole wavelength range, either resampled to 24 spectral channels or 
at the original spectral sampling. We then used two sets of PSFs to perform deconvolution of the cube at the different 
samplings. We use the deconvolution at the original sampling for extracting the spectra of the WR stars E2 and E4.0. 
The resampled deconvolution is used for E3.0 and 5.0. 

We also use the PSFs interpolated at the original sampling for creating a subtracted cube in which the four brightest 
stars (El, E2, E4.0 and the bright late-type star north-west of IRS13E) are subtracted. For this purpose we fit in 
each spectral plane the stars with Gaussians to determine their flux as a function of wavelength. Then we subtract 
the correspondingly scaled, interpolated PSFs from the original data cube. We use the resulting cube for analyzing 
the spectral features of the gas and the fainter stars. 

We flux-calibrate the SINFONI data by using the early-type star El, in the spectrum of which we manually replace 
the (few) spectral lines by continuum. The photometric errors from the SINFONI data are estimated by varying the 
selection of source pixels. During the fit of the such obtained data with an SED model, we rescale them such that the 
reduced = 1. 



Photometry from NACO images 

In order to measure the SED of the objects fainter or more confused than E5.0, we determine fluxes from the 
deconvolved broad-band and narrow-band images. We measure the flux of each object at each available wavelength by 
manually selecting the respective object pixels in the frames. For the broad-band images we estimate the flux errors 
from the standard deviation of the measured fluxes at the different epochs. For the narrow-band images we simply 
adopt the same flux errors as obtained from the broad-band images. 

For the broad-band images in H and Ks, we calibrate the magnitudes following iManess et al.l (|2007r FI. We derive 
the narrow-band calibration from the broad-band calibration by using the empirical extinction law 

= Aks * A/Aks" = Aks * A/Aks"'-'' (A1) 

from IPraind (|1989| ). As calibrators we use eight isolated, bright early-types stars with known temperatures within 
a distance of less than 2.5" to IRS13E. Since there is no suitable L'-band calibrator close to IRS13E, we use an 
indirect method by extrap olating the dereddened SED of the calibrators to L'-ban d and applyin g the extinction 
law from iLutz et al.l ()1996D longward of Brackett-/? and short-ward the power law of lDraini (Il989f). This results in 



AksM t,' = 1.75. The slope of the pow er law in equation I All is uncerta in: a = —1.6 (jRieke fc Lebofskvlll985t iRiekd 
Il999f) . a = -2 (jNishivama et al.ll2009D and a = -2 .21 fSchodel et al.l |2010D. Recent works may indicate that the 
extinction law used here is outdated |Stead fc Hoard [2009: Gosling ct al.l l2009f ). But at our precision the value of a 
does not affect the SED fitting because the shape of the derived SED is relatively insensitive to the value of a adopted. 
Nevertheless it introduces an uncertainty to the extinction values. Thus, it is important for the absolute luminosities 
of th e early-type stars. T o estimate these uncertainties, we use both the originally adopted extinctio n and the law 
from i Schodel et aTl (|2010l ) plus the magnitudes of the calibration stars in H, Ks and L' published by iSchodel et al.l 
(l2(m 

For the SED fits, we calculate the effective wavelength of the bands for each object. First, we apply an extinction 
value of Aks — 2.8 together with the atmospheric and the filter transmission to the source spectrum (either a blackbody 
of the output from the atmosphere modeling) . For the broad-band data it is necessary to perform a second iteration 
using the resulting SED in order to improve the effective wavelength. 

B: INDIVIDUAL OBJECTS IN IRS13E 
IRS13E1 

El is a hot OB supergiant (jPaumard et al.ir2006f ). Its fiux is constant in our data, and the spectrum does not reveal 
any signs of binarity. We use the H-band lines from the Brackett-series to determine the radial velocity (in K-band, 
the emission lines from nebular background are more dominant). We obtain fr.sp = 13 4±10km/s. This is 3a different 
to wlsr = 71 ± 20km/s, the Brackett-7 based value reported in iPaumard et all (|2006l ). 

We compare our atmosphere model of El to a simple blackbody model of T = 26000 ± 1000 K and do not see 
differences between the two models which would matter for the broad-band photometry. Therefore, we fi t the extinction 
using the photometric data and a blackbody model. We obtain Aks = 3.64±0.07 (figure[S|). Contrary to lMaillard et al.l 
(|200l . we do not find a weak dust component. This is likely because of the higher resolution of our data, since the 
SINFONI data show that the star is surrounded by gas probably also associated with dust. Fitting a supergiant 



' The calibration of lSchodel et al.l l|2010l ') is about 0.3 mag fainter 



in both bands. 
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Fig. 12. — Four of the 24 deconvolved narrowband images created from the high-quahty H+K-band SINFONI cube of IRS13E. 



model to the observed SED we obtain logL/L© — 6.23 ± 0.07. Thus El would be more luminous than all known 
supergiants of this spectral class. When using the calibration and extinction law of ISchodel et al.l (|2010D we obtain 
\ogL/LQ — 5.90 ± 0.07, consistent with the brightest known OB supergiants. 

IRS13E2 

E2 is a Wolf-Rayet star of type WN8 (jPaumard et al.l 120061 : iMartins et alll2007| ). There is no sign of binarity in its 
light curve or in its spectrum. The latter shows strong, broad emission lines. In addition, WR-stars have strong winds 
causing a near-infrared excess. Therefore it is necessary to use a full atmosphere model for the SED fit of E2. We 
manually adjust the temperature of a possible dust component and the extinction parameter to make the atmosphere 
model match the observed spectrum. 

The SED fit yields a second blackbody coniponen t of around T = 712 K and Aks = 3.51. This result is qualitatively 
consistent with the findings of iMaillard et alJ (|2004| ) . The smaller dust temperature they find is likely caused by their 
assumption that E2 can be modeled with a blackbody. We find that the dust component is significant in L'-b and, but 
in Ks-band the stellar component dominates (see Figure [3]), contrary to the proposal bv IMartins et al.l ()2007| ). 

We obtain for E2 logL/i© — 6.14 ± 0.1 when using the standard c alibration and e xtinct ion law of this paper and 
logL/L0 — 5.81 ± 0.1 when using calibration and extinction law of ISchodel et al.l ()2010[) . Thus. E2 is relatively 
luminous for its classification as WN8 star. E2 could be more massive and even younger (w 3 Myr) than most of the 
young, massive stars in the central parsec. However, due to the particular shape of the evolutiona r y trac ks, E2 could 
equally well have the same age as the other stars of similar type. In their figure 23 IMartins et al.l |2007ll gave for E2 
only an upper limit in luminosity because of the unknown dust contribution at the time. From our S ED fit we were 
able t o show that the Ks-luminosity is to about 96% of stellar origin. Therefore the luminosity of E2 in IMartins et al.l 
(|2007D is nearly identical with the value derived here. In any case, E2 is brighter than the other WN8 stars in the GC. 
This fact is is unaffected by calibration uncertainties. Also, binarity would not be able to fully explain the luminosity 
of E2. 

We use the 2.24 /xm NIII line from the full atmosphere model for deriving the radial velocity of E2 , which is 
more accurate compared to previous works. W e obtain ulsr = 65 ± 30km/s while IPaumard et al.l ((2006) reported 
'^LSR = 40 ± 40km/s and IMaillard et afl (|2004D gave wlsr — 30km/s. All these values are consistent with each other. 



IRS13E3.1 to IRS13E34 

The extracted PSFs of the SINFONI data are not good enough to measure the H-band part of the SEDs of the 
fainter objects around E3.0. Qualitatively, the spectra of E3.1 to E3.4 do not show any new spectral features compared 
to E3.0. There is only a gradual change in the spectral properties, such as the line velocities. Hence, these sources also 
likely consist mostly of gas and dust. This is confirmed by the SED fits, which are based on the narrow-band images. 
Because of the faintness of the sources we cannot reliably fit the extinction, which we restrict for pure dust sources to 
Aks = 3.4 - 3.9. 
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On most images in K-band (including the narrow-band images before 2005) and all L'-band data, the sources E3.1 
and E3.2 are strongly overlapping. Therefore we fit a combined SED for the two (figure [5]). The properties of the 
fitted blackbody are similar to E3.0 apart from its brightness. A stellar component is not necessary, the la limit is 
mxs =18.0, and thus again we can exclude that a WR star of the most common types is hidden in E3.1 and E3.2. We 
suspect that IMaillard et al.l (|200 4') found a stellar component in the eastern part of their source E3 due to the same 
bias as for E3.0 (section [4.2.ip . Also for the less well constrained objects E3.3 and E3.4 (table[2]) no stellar components 
are required. This conclusion does not depend on the assumption for the extinction. 

IRS13E3.5 

The SED of E3.5, derived photometrically, differs from the other sources in E3 (figure[5]). Its color is consistent with 
that of a star, and no dust is necessary to explain the SED of this object. Unfortunately, it is too confused and too 
faint for deriving its spectral type. Any temperature between 4000 K and 19000 K is possible, which also implies an 
uncertainty of A^ks = 0.4. The absolute magnitude of E3.5 corresponds to the red clump, which makes it likely that 
E3.5 is an old, late-type star. In addition, this star does not share the proper motion of the three brightest stars in 
IRS13E. We conclude that E3.5 most likely does not belong to IRS13E. 

IRS13E4.0 

E4.0 is a Wo l f-Ray et (WC9) star (jPaumard et a l."2006': 'Martins et al.l l2007D . Its radial velocity is very uncertain. 
IPaumard et al.l (j2006f) report ■^lsr = 56±70km/s, IMaillard et al. (2004) claim ■^lsr = ~30 km/s. From our atmosphere 
model (figure [3|) including wind induced line broadening we obtain vlsr — 200 ± 200 km/s. 

As for E2 we manually adjust the atmosphere model with an additional dust component to match the measured 
spectrum (figure[3]). Clearly two components are needed to describe the SED, but the derived values are very uncertain. 
The dust component with T = 1400 K is significant. It contributes 38% of the flux in the Ks-ba nd, as indicated by 
iMartins et al.l (|2007[ ). Our measurement is more precise than the one from IMaillard et all ()2004D . because we use a 
full atmosphere model and extend the wavelength coverage to L'-band. We obtain the same extinction as for El. The 
luminosity of E4.0 is logL/i© = 5.23. Because of the difficult modeling of this star and WC9 stars in general it is not 
possible to conclude whether E4.0 has an age similar to the other WC9 in the GC. 

IRS13E4.1 to IRS13E4.3 

Besides H- and Ks-band, E4.1 is detected also in L'-band and has overall a red SED. It is too close to E4.0 to 
allow the broadband SED to be extracted from the SINFONI data. The spectrum shows only gas lines and no stellar 
features. For fitting the SED it was necessary to fix the extinction at the upper end of the possible range to ^ks = 3.9. 
With this the data are fit well by a single blackbody of T = 844 ± 24 K. Any stellar component has to be fainter than 
TTlKs = 18.4. 

E4.2 does not show any spectral features and is not detected in L'-band. Therefore, there is less data than for the 
previous sources. We fix Aks = 3.9. From a single blackbody fit we get a temperature of less than 2800 K, which is too 
cold for a star of mxs = 15.4. Thus, we need to add a second blackbody to the fit, which means we classify E4.2 as a 
star. As in the case for E3.5, this faint star matches the typical luminosity of red clump stars. Hence, it is most likely 
an old late-type star. Its motion, being slower westward than for the three bright stars in IRS13E, argues against E4.2 
being a physical member of IRS13E. 

For E4.3 we are only able to photometry in Ks-band, yielding mxs — 15.5 ± 0.1. 

IRS13E5.0 

The source E5.0 is isolated enough to extract its broadband SED in the Ks-band from the SINFONI data (figure [S]), 
but because of its location close to the edge of the data cube, the SINFONI H-band measurements are not reliable. 
Hence, we combine the NACO-based fluxes with the SINFONI data. We allow for a global scaling factor to match the 
two data sets, which yields an offset of 0.24 mag. In any case, the flt does not change significantly by either including or 
excluding the SINFONI data. The preferred fit yields Aks = 3.22 and a dusty blackbody. Thus, we fix the extinction 
value to the lowest value from the range assumed here to Aks = 3.4. Again, no stellar component is needed, and any 
star inside i s fainter than mxs = 18.5. The dust of E5.0 is with T = 700 ± 8K colder than for most other objects in 
E3 and E4. IMaillard et al.l ()2004D were only able to identify one source in E5. The scatter in their data points is so 
high that it is not clear how significant the stellar component is in their fits. 

IRS13E5.1 to IRS13E5.4 

The other sources of E5 are also red but fainter than E5.0. They appear embedded in the brightest part of the 
minispiral. It is not possible to measure their SEDs separated from the minispiral at our resolution, because the 
background is inhomogeneous and therefore the background subtraction too uncertain. This holds both for the NACO 
photometry and for the SINFONI data. In addition, the sources E5.2 and E5.3 are located on the deconvolution ring 
of the bright source El. Thus, the fluxes of these two sources are biased to high values, especially in H-band. Finally 
many of the objects are not detected in all our NACO bands. Consequently, the characterization of these source is 
considerably more uncertain than for the other sources. The spectra of these sources show only gaseous emission 
lines. We again assume that the continuum emission is mainly caused by dust, and thus we restrict Aks to the range 
3.4-3.9. 
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Fig. 13.— SINFONI H+K-band spectrum of E6. 



The source E5.1 is a double source on the higher quahty Ks-band images obtained since 2008, but we are not yet 
able to disentangle them. This source is also detected in the L'-band. With our assumption for the extinction, the 
best fit is a blackbody with T = 947 ± 25 K; any stellar component needs to be fainter than ttiks = 17.8. 

The source E5.2 is not detected in H-band but in L'-band, where it is nearly overlapping with E5.3. We nevertheless 
simply assume that all of the L'-flux is caused by E5.2. With this we obtain a temperature of T = 642 ± 70K from 
the SED fit, which is consistent with the non-detection in the H-band. Any star hidden in E5.2 must be fainter than 
TTiKs = 15.6. Given the lower temperature, E5.2 is likely a part of the minispiral. 

Contrary to E5.2, the object E5.3 is not detected in L'-band but in H-band. Our fit does not require a stellar 
component. The Icr-limit is toks = 16.0 and we obtain a temperature of T = 978 ± 524 K. The large error is caused 
by the smaller wavelength range on which the SED fit is based. 

The source E5.4 is located close to E5.0. It is detected only in Ks- and the narrow bands. Hence we fit only a single 
blackbody with the restricted extinction interval. We obtain a temperature of T = 817 ± 474 K. 

IRS13E6 

The spectrum f figure [1^)1 of E6 {niKs = 13.82 ±0.10) shows CO-ba nd heads in H- and K s-band and lines of NaT and 
Cal. Thus, E6 is a late-type star with a spectral class around K3HI. iMaillard et al.l (|2004[) identified E6 a s an -star, 
but they used data from three relatively broad bands for their photometric identification. iBuchholz et al.l ()2009f ) used 
narrow bands around the CO-band heads and determined correctly that E6 is a late-type star. We derive a radial 
velocity of wlsr = 106 ± lOkm/s from the CO-band heads in Ks-band. 

In order to avoid biases, we exclude from the SED fit the bands which contain flux of the CO band heads. The fitted 
extinction value of — 3.68 ± 0.09 is consistent with the extinction of El. Due to its faintness E6 cannot be a red 
supergiant and it is thus older than the three bright early-type stars. 



C: ASTROMETRIC DATA 
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Fig. 14. — Positions and linear motion fits for the objects in IRS13E. Top: The three bright stars in Ks-band (asterisks and sohd fines) 
and H-band (boxes and dotted fines). Bottom: Some of the fainter objects. All data and fits are offset from each other for clarity. 



